Correlations in non-equilibrium Luttinger liquid and singular Fredholm determinants 
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We study interaction-induced correlations in Luttinger liquid with multiple Fermi edges. Many- 
particle correlation functions are expressed in terms of Fredholm determinants det(l + AB), where 
A(e) and B(t) have multiple discontinuities in energy and time spaces. Such determinants are a 
generalization of Toeplitz determinants with Fisher-Hartwig singularities. We propose a general 
asymptotic formula for this class of determinants and provide analytical and numerical support to 
this conjecture. This allows us to establish non-equilibrium power-law singularities of many-particle 
correlation functions. As an example, we calculate a two-particle distribution function characterizing 
correlations between left- and right-moving fermions that have left the interaction region. 

PACS numbers: 73.23.-b, 73.50-Td 
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Non-equilibrium phenomena in (effectively) one- 
dimensional correlated systems — including Kondo and 
related quantum impurity models 1 - - — , quantum Hall 
edge interferometry&~— and energy relaxatio n 13 ! 14 , car- 
bon nanotube tunneling spectroscopy^ 5 -, Fermi edge 
singularity 1 -, and correlated electrons in quantum 
wires^"— — are attracting lots of experimental and theo- 
retical interest. In these problems applied voltages lead 
to formation of distribution functions with two or more 
Fermi edges, which generates non-equilibrium scaling, 
criticality and decoherence. A remarkable property of 
the Luttinger liquid (LL) model is a possibility of ex- 
act solution even for such non-equilibrium distributions. 
In this paper we show how many-particle fermionic cor- 



relation functions G r . 



(*}... ... *2n) can be 

calculated and discuss the underlying physics. It turns 
out that G n are given by a certain type of singular Fred- 
holm determinants. Our result for the asymptotics of 
such determinants is expected to be relevant also to other 
non-equilibrium many-body problems. 

In previous works 1 ^— , we have shown that single- 
particle Green function (GF) in the case of LL (and in 
a number of related problems) can be expressed through 
Fredholm determinants of a "counting" operator 



A [A, B] =det(l + AB), 

-iS(t) 



(1) 



where operators A(t) = e w— 1 and B(e) — n(e) are 
diagonal in the time and energy representations, respec- 
tively, and [e, t\ = i. The one-particle distribution func- 
tion n(e) characterizes the non-equilibrium state of in- 
coming electrons, while the time-dependent phase S(t) 
encodes the information about the interaction. For a LL 
adiabaticaly connected to reservoirs the phase is 1 ^ 



6(t) = S 1 e(\T\/2-t) signr, 



(2) 



where 0(t) is the Heaviside step function. This allows 
one to reduce the operator in Eq. ((T|) to a Toeplitz 



matrbsSi, with 1 + n(e)(e~ lSl — 1) determining its sym- 
bol. (For a non-adiabatic coupling one gets a sequence 
of rectangular pulses in 6(t), yielding in the long- wire 
limit a product of Toeplitz determinants.) When n(e) 
has discontinuities ("Fermi edges"), the single particle 
GF acquires a non-trivial power-law behavior. This is in 
particular the case for multi-step distributions 



n(e) 



1 = a , 
ai , 

dm ) 

= a-m+l , 



e < e 

e < e < ei 



(3) 



_i < e < e r) 
< e . 



The low-energy behavior of the single-particle correlation 
functions can be understood with the help of the general- 
ized versio n 21 ' 23 of the Fisher-Hartwig conjecture^ (see 
also 25 ' 26 ). Under generic conditions, it yields a power- 
law energy dependence masked by dephasin g 21 ' 23 . When 
the electronic system is in a pure state, i.e. for ([3]) with 
all a, equal or 1, the dephasing is absent and the effects 
of correlations are particularly strongly pronounced. 

Higher-order GF of a non-equilibrium LL can be cast in 
the form Eq.([T]) as well, which offers a remarkable oppor- 
tunity of exploring many-particle correlations in a quan- 
tum many-body system out of equilibrium. However, 
there is a serious complication as compared to the single- 
particle GF case. Since each creation or annihilation of 
an electron induces a jump in S(t), the latter has a form 



S(t) 



= <5 , 
or, 

o = 4+i 



t < t Q 

t < t < tl 

tk-1 <t <tk 

tk <t. 



(4) 



As we see, the determinant is now generated by func- 
tions A(t) and B{e) that both have multiple jumps and 
therefore is not of Toeplitz type. Since time and energy 
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FIG. 1: Construction of the matrix fiij, Eq. ©, determining 
the asymptotic behavior of the Fredholm determinant. 



enter Eq. ([TJ on equal footing, one may expect the t and 
e discontinuities to play a similar role in the asymptotic 
behaviour. This suggests that there should exist a gen- 
eralization of the Fisher-Hartwig formula valid for this 
class of determinants. This generalization represents one 
of key results of the present paper. 

We first state our result for the long-time behavior of 
the determinant ([1} and then present arguments in favor 
of it. To formulate the result, it is convenient to draw 
a grid in the time-energy plane defined by points where 
A(t) and B(e) exhibit jumps, as shown in Fig. [T] (for the 
case k = 2,m = 2). This divides the plane in (k + 2) x 
(to + 2) rectangles (with the "external" ones extending 
to infinity) . With each of them we associate a number 

La = (i/27r) log [1 + (e-* - l)a 3 ] ; 

< z < fc + 1 , 0<j'<to+1. (5) 

We impose the condition that a branch of the logarithm 
at infinity is fixed by L^ = 6i/2ir and L^ m+ i = Lqj = 
Lfc+ij = 0, while for finite pieces any branch can be 
chosen. We now define a matrix 



(6) 



where indices i = 0, . . . , k and j = 0, . . . m correspond to 
the grid lines 2 ^, and a set of exponents 



— Pij Pi' j' Pi]' j v , 



(7) 



The asymptotic behaviour of the normalized (to its 
zero-temperature form) determinant A[5(t), n(e)] = 



A[6(t),n(e)]/A[6(t),T = 0] is given by 
A[S(t),n(e)} = ^C[ y 8]exp(-i^t i ft j e i ) 



br. 



x n n [fe-Mfe-^)p'^.(8) 

i.i' — j,j'=0 
i'<i j'<j 

Here the sum 5^ br runs over all possible branches of the 
logarithms in (|5|), and C[j3] are numerical coefficients 28 . 
Equation ([8]) represents the infrared asymptotics which 
holds provided that (ti — t -)(ej — e'j) ^> 1 for all i > i' and 
j > j'. If for some set (i,i',j,j') the opposite inequality 
holds, the corresponding factor should be dropped in ([5]). 
When all inequalities (ti — t v ) (ej — ey ) ^> 1 are fulfilled, 
one can calculate the total power of each of the factors 
ti—tii and ej — eji in (|8j) by using the sum rules J2i Pij = 
and £\ P^ = (5 i+1 - 6^/2%, which yields 



A = ^C[/3]exp[- ? 5>/%(6,-A)] 

br. i ,j 

x Y[[A(t^t,)r" UKej-ej^/Api', (9) 

i'<i j'<j 

where [i lV = J2±PijPi'ji Vjj> = J2iP*jPij'i and A is tne 
ultraviolet cutoff. 

While we have no mathematical proof of Eq. (|5J|, we 
have strong evidence in favor of its validity. Heuristically, 
Eq. ijSJ) represents a natural extension of the generalized 
Fisher-Hartwig formula of Refsi 21 ' 23 (valid for Toeplitz 
determinants) onto the present case. Indeed, the Fisher- 
Hartwig formula suggests that power-law factors in the 
asymptotics of Toeplitz determinants are due to the pres- 
ence of singular points (ti, Cj) and each pair of such points 
contributes independently to the result. Physically, the 
contribution of the pair of points (ti, e,-) and (ty, 6j>) rep- 
resents the effect of particle-hole pair scattering from the 
vicinity of the Fermi edge €j to Fermi edge ej/ within the 
time ti — tii. Combined with the symmetry between e 
and t, these arguments naturally lead to Eq. ©. On a 
technical level, we have proven a generalization of strong 
Szego limit theore m 29 ' 30 to the case of a multistep A(t) 
and smooth B(e). Substituting in this formula a multi- 
step -B(e), we obtain logarithmic divergencies in the ex- 
ponent, yielding power laws of Eq. ©. We refer the 
reader to Supplementary Material 32 for detail. Finally, 
we evaluated the determinant numerically for the sim- 
plest non- Toeplitz case, with 8(t) having three jumps. 
The results^ 2 - confirm Eq. ([5]). 

We now apply Eq.® to the problem of correlations 
in a non-equilibrium LL. We consider a LL conductor 
of length I (characterized by LL parameter K) driven 
out of equilibrium via the injection of electrons with dis- 
tributions n^ n (e) [rj = R (L) for right- (left-) movers] 
from the non-interacting leads, see Fig. 1 of Refj^. To 
be specific, we assume that the right-moving electrons 
have the triple-step distribution function with popula- 
tion inversion (Fig. [5]) while the left movers are at zero 
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FIG. 2: Distribution functions of the incoming (top) and out- 
going (bottom) electrons. The incoming right-mover distri- 
bution function has three Fermi edges and exhibits inverse 
population. Incoming left movers are assumed to have zero- 
temperature Fermi-Dirac distribution. Distribution functions 
of fermions at the output of the device show power-law sin- 
gularities at Fermi edges with exponents indicated in the leg- 
ends. To generate the plot we assumed K — 0.8 
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FIG. 3: Top and Middle. Phases Srr (£) and &RL(t) control- 
ling the distributions of left- and right- mover at the out- 
put of the wire. The time interval 2KI/vf between dif- 
ferent pulses is macroscopically large and the correspond- 
ing determinants can be approximated as products. Bottom. 
The phase 5n(t) — Srr (£) + 5rl(£) controlling the correlator 
{{riR^XR, tR)nh{xL, el))) at xr + xl = —Kl. 



temperature. This is the simplest non-equilibrium setup 
without dephasing. The width e 2 — eo of n x R {e) is set 
to unity, thus determining the time and energy scale 
of the problem. We model the leads by ID conductors 
with K — 0. The correlations effects discussed in this 
paper can be traced back to the scattering of LL plas- 
mons at the boundaries of the wire. This scattering is 
strong when boudaries are sharp (compared to the plas- 
mon wave length), K(x) = KQ(l/2 — \x\), so that we 
focus on this regime. We stress that we assume the ab- 
sence of the electron backscattering in the wire, i.e., K{x) 
is smooth on the scale of the Fermi wavelength. 

At thermal equilibrium (n 1 ^ L (e) are Fermi-Dirac dis- 
tributions with equal temperatures), the interaction 
causes no correlations for electrons outside the LL wire. 
In other words, electrons leave the wire in the same state 
as free fermions would do at given temeparture. The 



situation changes dramatically under non-equilibrium 
conditions^!. The plasmon scattering at the bound- 
aries of the wire (characterized by reflection coefficient 
TZ = (1 — K)/(l + K) and transmission coefficient T = 
y/l — TV 2 ) not only leads to an electron energy redistri- 
bution but also induces correlations between outgoing 
electrons. As we found earlier—, distribution functions 
n< RL( T ) m the non-interacting regions are given (up to 
a numerical factor) by the determinant (TT]) with phases 
S RR (t) and o~L R (t) shown in Fig. [3] Both phases consist 
of an infinite sequence of rectangular pulses separated 
by a time interval 2KI/vf- In the long- wire limit the 
corresponding determinant can be decomposed into the 
product of Toeplitz determinants controlled by individual 
pulses 3 - 3 -. An asymptotic analysis of these determinants 23 
leads to the conclusion that distribution functions of out- 
going electrons have power-law singularities at the Fermi 
edges (Fig. [2]). Note that for relatively weak interac- 
tion the interaction effects are most pronounced at the 
inverted edge of the n R (e) signaling that the dominant 
physical process is the scattering of electron-hole pairs 
from this edge to the Fermi edge of ul (e) . 

To reveal the correlations induced by the plasmon scat- 
tering, we consider the irreducible two-particle distri- 
bution function of outgoing electrons in the two leads, 
f(e R ,e L ,x R +x L ) = ((n R (xR,e R )h L (x L ,e L ))) with x R > 
1/2 and xl < —1/2. It is a non-trivial function of x r +xl, 
with the correlations being most pronounced in a vicinity 
of the points x R + xl = (2m + l)Kl 3 K We focus on the 
case x r + xl = —Kl where the correlations are maximal. 
At this point the function f{e R ,eL,—Kl) = f—\(e R ,eL) 
is given by the determinant ([1]) with the phase S R (t) = 
S RR (t) + 5 LR {t), see Fig. H Explicitly 



f-i(eR,e L ) 



1 



dr R dT L 



,i<iRTR+ie L T L 



47T 2 J T R T L 

(K[6 R , n£] - A[S RR , n}R]A[S RL , n£]).(10) 



Investigation of the correlations of left- and right- movers 
is now reduced to calculation of the Fredholm deter- 
minant with phase 6 R (t) that can be readily done by 
employing Eq. ([5]). For a long wire the determinant 
A[<5/{, ri 1 ^] decomposes into a product of determinants 
corresponding to individual pulses forming S R (t) (see 
Fig. [3]). We write A[5 R ,n R ^} = A{A 2 , where Ai is a 
contribution of the first pulse and A2 of the remaining 
pulses. For definiteness, we focus on the case of a weak 
interaction (small TV) when the phase 5 R (t) is small out- 
side the first pulse. Accordingly, the asymptotic behavior 
of the determinants for all pulses but the first one is gov- 
erned by a single term in the sum (j8]) corresponding to a 
choice of logarithm branches such that | L%j \ < tt for all i 
and j. [The analysis for an arbitrary interaction proceeds 
in the same way; one just may need to take into account 
several further terms in Eq. ([8]).] This yields 



A 2 = C 



( TR+T L 
V T R - TL 



(11) 



4 




FIG. 4: Left and Middle. Contour plot of the irreducible correlation function /-i(tr, tl) multiplied by trtl- The left panel is 
the result of numerical evaluation of Fredholm determinants. The middle panel shows the fit based on Egs. (|ll[l and (1121) with 
coefficients d used as fitting parameters. Shadowed regions on both graphs are outside the applicability of the asymptotic 
expansion ([5|. Right. Irreducible correlator of fermionic occupation numbers in energy representation. 



To find Ai, we note that SR,(t) is close to 2tt; an in- 
spection of Eq. © leads to the following three dominant 
contributions: 



Ai = de- 



Co e" 



f T R + T L 



X T 



-2n 2 (K 2 +2) _ 2U 2 



(12) 



Since we are not interested in the dependence on ener- 
gies ti — €j, we have absorbed it in the coefficients Cj in 
Eqs. (11]) and (Jl2j). Combining (1111) and (1121) with asymp- 
totic expansion for the Toeplitz determinants A[(5r^] and 
A[<5al}^, we arrive at the asymptotic expansion of irre- 
ducible correlation function /_i in time domain. 

The two-particle correlation function /_i(ta, tl) is 
shown in Fig. 2] for K — 0.8. The left panel presents 
the results of direct numerical evaluation of the Fred- 
holm determinants, while the middle panel displays the 
asymptotics given by Eqs. (fTTj) and (fT2|) with the coef- 
ficients Ci used as fitting parameters. As expected, the 
asymptotics based on Eq. (|8]) reproduces correctly the be- 
havior of /_x(tr, tl) outside the region of (tr, T£)-plane 
where tr < 1 or tl < 1, or \tr — tl\ < 1. 

Power-law long-time behavior translates into singular- 
ities in the energy reperentation of the correlation func- 
tion /_i(e,R, £l) plotted in right panel of Figj4] Strong 
correlations are observed that form a "crest" going along 
the line 6r + €l = ei (thick red line). The analy- 
sis shows that the main contribution to this structure 
in /_i(efl,£i) comes from the vicinity of the main di- 



agonal \tr — tl\ < 1 of the (tr, T£)-plane. In this 
region, the main contribution to Eq. (fl2"|) is given by 



the last term yielding Ai 



'(tr + T L ) ai with 



oj = 4ft - 6ft 2 + 4ft 3 -2ft 4 , whereas Eq. reduces to 
A 2 ~ {tr + tlT 2 with a 2 = -2T 2 ft 2 (l-ft) 2 /(l+^ 2 )- 
This yields the behavior of f-i(eR, €l) near the crest 



(13) 



Physically, the crest originates from intensive scattering 
of electron-hole pairs from the vicinity of the left-mover 
Fermi surface e = to the inverted Fermi surface ei of 
right movers, and Eq. ([TB")) can be viewed as a Fermi-edge 
singularity in a higher-order correlation function. 

To summarize, we have studied correlation induced by 
interaction in a non-equilibrium LL. We have shown that 
electron scattering between multiple Fermi edges leads 
to power-law singularities in many-particle distribution 
functions. As a particular example, we calculated a two- 
particle distribution function characterizing correlations 
between left- and right-moving outgoing fermions. Tech- 
nically, many-particle correlation functions are expressed 
in terms of Fredholm determinants det(l + AB), where 
A(e) and B(t) have multiple discontinuities in energy and 
time spaces, respectively. We have conjectured a gen- 
eral asymptotic formula for this class of determinants and 
provided an ample analytical and numerical support to 
this conjecture. The results are expected to be relevant 
to a broad class of non-equilibrium many-body problems 
including the Kondo problem; work in this direction is 
currently underway. 

We acknowledge support by Alexander von Humboldt 
Foundation, ISF, and GIF. 
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Appendix A: Analytical justification of the asymptotic formula for singular Fredholm determinants 

In this section we present analytical arguments in favor of the Eqs. (8), (9) of the main text for the asymptotic 
behavior of the determinants det(l + AB) with A(e) and B(t) having multiple discontinuities in energy and time 
spaces, respectively. We begin [Sec. IA 1| by reminding the reader about Szego formula for Toeplitz matrices with 
smooth symbols and a generalized version of the Fisher-Hartwig formula for Toeplitz matrices with singular symbols. 
We emphasize there a connection between these two formulas. We proceed then by formulating and presenting a proof 
of a generalization of Szego formula for the case of a multi-step function A(t) and smooth B(e), Sec. IA 21 Finally, in 
Sec. IA 31 we use this result to obtain the scaling behavior of Fredholm determinants in the case when both functions 
A(t) and B(e) have multiple singularitites. 



1. Toepliz determinants: Strong Szego limit theorem and generalized Fisher-Hartwig conjecture 

Let us consider the determinant of N x N Toeplitz matrix Tjv[/] = {.fi-j}, < i, j < N — 1 generated by function 
f(z) defined on the unit circle z = e and having sufficiently smooth logarithm V(z) = ln/(z). We use the notation 

/W=//W^-^. (AD 
Smoothness of V(z) on the unit circle implies the existence of the Wiener-Hopf decomposition for / 

f(z) = e v + (*) e Vo e V-(*) , Vo = 1 l m *± (A2) 

Zm J z 

with functions V±{z) being analytic inside and outside the unit circle respectively. In terms of Fourier components of 
V(z) 

V + (z) = jrV k z k , V-(z) = jrV- k z- k . (A3) 

k=l k=l 

The strong Szego limit theorem— states that under the assumptions made above the asymptotic behavior of the 
Toeplitz determinant det T/v [/] at N — > oo is given by 



detTjv[/] = e NVo exp 



Y, k VkV- k 



Lfc=l 



■o(l) . (A4) 



When interpreted in physical terms Eq. (|A4[) yields the long-time asymptotic of the single-particle correlation functions 
of II? interacting fermions in a non-equilibrium state characterized by smooth distribution function^. The exponential 
term in (|A4[) encodes herewith the information about the oscillations of the Green functions and its exponential decay 
due to non-equilibrium dephasing, while the precise value of the pre-exponential factor is of little importance. 

The situation changes significantly if one considers single-particle Green functions in a state characterized by 
distribution function of electrons having m + 1 = 1,2,... discontinuities (Fermi edges). To be specific we limit our 
discussion to the case of multi-step distribution function 

1 = a , e < e 

ai , e < e < £i 

n«)={ ... (A5) 

= a m+ i , e m < e , 

Supplementing the theory with the ultraviolet cutoff A, one ends up with a Toeplitz matrix with symbol having 
m + 1 jumps. In the context of Toeplitz matrices such discontinuitites are known as Fisher-Hartwig singularities. 
More specifically, the symbol of interest is 

rn 

f(z) = e v °z^T=o ft Y[ g Zj 0j ( z)z "ft , (A6) 
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FIG. 5: Pictorial representation of Eq. (|A13|l . A set of numbers fin is associated with the set of crossing points (ei,ti) of 
singularities in (e, t) plane. Each pair of such crossing points with different time coordinates gives rise to a logarithmically 
diverging contribution to the exponent in Eq. (|A13[) . 



with Zj = e lnej / A = e l6j and 



9*i ,ft( z ) = i P -»rft- 



-7T < argz < 0j 
6j < are z < ir . 



Here the numbers /3j controlling the jumps of V(z) = ]nf(z) at points Zj are given by 



1 



In- 



2ni 1 



(A8) 



(A9) 



We are interested in the large- TV asymptotic behavior of the determinant det Tn [/] . In physical language the size of 
the Toeplitz matrix corresponds to the time in the correlation function r = ttN / A. 

Strictly speaking, the strong Szego theorem does not apply to the Toeplitz matrix with a singular symbol (|A6[) ; 
the corresponding extension of the large- N limit theorem is known as Fisher-Hartwig conjecture^. This conjecture 
states that at large N 



detT N [f] = C[p]e NVo N-ZT=o% J] \ Zi _ z . 

0<j<£<m 



2ft ft 



(A10) 



with C[/3] being a (known) numerical coefficient. Casting this result in terms of physical variables and assuming that 
all energies q are small compared to A, one obtains 



detT N [f] =C[P} ( 



L +*5ft"L ft 



(JL 



77 \ Ift"Lo Pj 



n g(. 



(All) 



0<j<i<m 



While a rigorous mathematical proof of Eq. (|All[ ) is highly non-trivial and for the general settings was achieved only 
very recently, a strong argument in favor of Eq. (|Alip can be given on the basis of the Szego formula (|A4[) . Indeed, the 
jumps of V(z) at points z, = e l6i dictate the following asymptotic behavior of the Fourier coefficients Vj. at k —} oo 



i2> 



(A12) 



The naive application of the Szego formula leads now to 



detTjv[/] = e w °exp 



m oo ^ 



»J=0 



fe=l 



(A13) 



where we have introduced here (3u 
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Equation (|A13I) can be given a pictorial representation shown in Figj5] We attribute the quantities 0u to the points 
(ej,i() in (e,t)-plane where jumps of the functions n(e) and 5(t) cross. Each pair of such points with different time 
arguments [i.e., {e^ti) and (a',ti>) for t\ ^ ty)] provides the following contribution to the last factor in Eq. (|A13|) : 



exp 



>uPi> 



E k 

k=l 











~ exp 









(A14) 



To give meaning to this expression, we have to cut off the logarithmically diverging integral. The lower integration 
limit is set naturally by the ultraviolet cutoff A. For i — i' one expect the upper integration limit to be given by 
the time r. On the other hand, for i ^ i' the upper limit of integration is expected to be l/|ei — (recall that the 
asymptotics we are discussing is valid under condition r|ej — | 3> 1). Using these simple cutoff rules together with 
the expression (IA7|) for the coefficient Vb, we see that the Szego theorem justifies in a natural way the Fisher- Hartwig 
conjecture (|A11I) . 

There is, however, the following subtlety here. Recalling the definition (| A6|) of the symbol with Fisher-Hartwig 



singularities, one sees that the numbers f3j are defined modulo a set of integer numbers n 3 - such that 



E 1 

3=0 



= 0. 



(A15) 



In other words two sets of coefficients {/3?} and {f3j = /3j + rij) describe the same symbol f(z) and should produce 
the same asymptotics. The recently proven rigorous formulation of the Fisher-Hartwig conjecture^ indeed respects 
this requirement. Specifically, it states that one should choose the set of f3j in Eq. (| Al 1|) in such a way that the sum 
S j=o Pj determining the power of N takes the smallest possible value. 

More recently, it was shown i n 21 ' 23 that other choices of integer shifts (i.e. of logarithm branches) in the expressions 
for Pj are also important and determine the sub-leading terms in the asymptotic behavior of a Toeplitz determinant 
with a singular symbol. A generalized Fisher-Hartwig conjecture formulated in this work reads 



dctT Ar [/]=^C[/3]e^ 
{0} 



n 

0<i<i<m 



2ft ft 



(1 + ...). 



(A16) 



The summation here goes over all choices of {/3j}, i.e., over all sets of integer shifts nj satisfying Eq. (IA15|) . Each such 
set yields in Eq. (|A16|) a term with a distinct oscillatory exponent e l ^i= Pi e 3 T , Equation (IA16P presents explicitly the 
leading asymptotic behavior for the factor multiplying each of these exponents. Apart from this dominant term, there 
will be in general also subleading (in powers of 1/t) terms corresponding to the same exponent; these are abbreviated 
by + ... in the last bracket. 

While a rigorous proof of the generalized Fisher-Hartwig conjecture Eq. (|A16I) remains to be constructed, there is 
little doubt that it is correct. 



2. Multiple singularities in A(t): Generalization of Szego formula 

We assume now that the electronic distribution n(e) is a smooth function of energy, whereas the scattering phases 
S(t) has a multi-step structure: 



5(t) = 



= S , 
Si, 

Sk , 

o = 4+i 



t < t 

t < t < t 1 

tk~l <t <tf- 
t k <t. 



(A17) 



We will assume for definiteness that k = 2; generalization to a larger number of steps in S(t) is completely straight- 
forward. Among known derivations of Szego formula, the operator-theory approach due to Widom 3 ^ permits a 
particularly transparent generalization to our problem. We thus follow this approach here. 
Upon introduction of an ultraviolet cutoff A the functional determinant of interest 



det 



1+ e 



(e-^ - l) n(e) 



(A18) 



is reduced to the determinant of a finite matrix of the size N x N 



Tni,N 

v Tn 2 ,n 



"/(i) 



with the structure 
, N = N 1 +N 2 . 



■ Ni x N with a symbol 



Here Tjv 4) jv [/^l stands for a rectangular Toeplitz matrix of the size 

= [1 + (e-** - 1) n(e)] e "^/ 2A 
which is smooth on the unit circle z — e lTrc / A . The numbers Nq and N\ are determined by S(t) via 

iVi = (ti-to)A/7r, N 2 = (t 2 -t 1 )A/ir. 
Let us introduce a semi-infinite Hankel matrix with symbol / according to 

H lf] = {fi+j+i} > i, j = 0, . . . oo 
and the semi-infinite matrices Pn and Qn with matrix elements (i, j — 0, . . . oo) 

max(i, j) < N — 1 



, \ _(Sij, max(i,j) < N- 1 , fn > _jk,N-i- 1 , max(i, j) < N - 1 
W« ~ \ 0, max(i,i) > JV- 1 a Q WArj y ~ \ , max(i, j) > JV — 1 . 

With this notations matrix An ,Ni [/ i f^] can be presented in the form 



An!,n 2 



,(2)1 = ( t Ni [/«] Q^iT [/( . 
1 " "[f (2) ]Q Nl T N2 [frt] 



Here the symbol / is obtained from a symbol / according to 

f(z) = f(z- 1 ) . 

The block representation (|A24[) of the matrix An 0i n 1 allows us to write its determinant in the form 
detA NuN2 = dct [T Nl [/«]] det [t N2 [f^J\ det f - T Na [f^ ^ H [f^ T Nl [/«] ~* H [/«] P N , 



(A19) 

(A20) 
(A21) 
(A22) 

(A23) 

(A24) 

(A25) 
(A26) 



We have used here that QnTn[/]Qn = It is important that the last determinant in Eg. IA26l has a finite limit 

as Nq and N\ go to infinity, 



C = dct 



1 — T 



(2) 



H 



(2) 



T 



U) 



H 



(A27) 



Here T[f] stands for semi- infinite Toeplitz matrix with symbol /. We can now use simple properties of semi-infinite 
Toeplitz and Hankel matrices^ 



H[ < t>]T[i,]+T[ct>]H[^] = H[^\ 
to bring the expression for constant C to the form 



C = det 



T 



(2) 



T 



/ (2) // (1) l T \l/fW 



On the other hand, one can recast the Szego result (|A4I) into the forrn^ 

detT N [f] = e NV °det [Ttf-^Ttf]] . 
Combining now Eqs. (|A30[) and (|A31[) we get for the determinant of interest 

det = eW^+W™ x det W \l/.f^] T \.f {2) / f {1) ] T , 



(A28) 
(A29) 



(A30) 



(A31) 



(A32) 
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where Vq is the zero Fourier component of /W. The derivation above can be immediately generalized to the case of 
S(t) having arbitrary number of steps, Eq. (|A17l) . The resulting determinant reads 



detA 



N,,...N k 



with 



/(i+D 



3 Ef =1 W D w xdet 



T 



,(*0 



.T 



,(o) 



.9 



i = 0, . . . fc and 

f(p) = j{k+i) = x _ 



Note that by definition 



(A33) 

(A34) 
(A35) 

(A36) 



i=l 



Equation (|A33[) allows to reduce the problem of asymptotic behavior of det AN lt ...N k to the evaluation of a deter- 
minant of a product of semi-infinite Toeplitz matrices. This is a crucial step forward because it enables us to apply 
the Winer-Hopf method to the problem. Specifically, let us introduce the Wiener- Hopf decomposition of the functions 



7 (0 



gV(z) = e v ?We u ie u -M 
e unit circle. The analytic 

T[ 9 W]=e< ) e T { U -]e T l U +] 



Here Utys (z) are analytic inside (outside) the unit circle. The analytic properties of U± allow us to write 



Thus 



det .4 



Ni,...N k 



(fe) 



= e £iU"«vtf« xdet 



M, k) e T [ U - k) ] e T [K^] P U 1 ^ T [u™ } T [U* 



Note that due to the relation (|A36[) 



K K 

E^ = E r [ 



u 



= 0. 



(A37) 



(A38) 



(A39) 



(A40) 



and one could naively conclude that the determinant in (|A39I) is equal to unity. This is not true, however, due to the 
infinite size of the matrices involved. Instead, one should make use of the formula^ 



det e UR ... e Ul = 



exp 



l<j<i<R 



(A41) 



valid for arbitrary set of operators Ui satisfying X)iLi ^ = 0- Here [*, *] stands for the commutator. Applying (IA41[) 
to (|A39[) . we obtain 



det ,4 



Ni,...N k 



(fe) 



x exp 



E tr 

0<i<j<k 



T 



U 



U) 



■ T 



U 



(i) 



(A42) 



Finally, evaluating the traces one finds 



det A 



Ni,...N k 



f {1 \-.-f 



(fe) 



x exp 



- E J2^ j)u -l 

0<i<j<k q=l 



(A43) 



Equation (|A43|) is the generalization of the Szego limit theorem to the case of determinants with a multi-step A(t) 
[and smooth B{e)\ and constitutes the main result of this subsection. In the next subsection we will analyze its 
implications for the case when B(e) is not smooth but rather has multiple steps. 
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FIG. 6: Construction of the numbers fii v determining the asymptotic behavior of the Fredholm determinant. Each of the 
coefficients /3i p can be represented as /3ij = Li.j + Li+i.j+i — — with Ly = ;p log [l + (e~ lS * — l)a 3 -l . 

3. Multiple singularities in A(t) and in B(e): Further generalization of Fisher-Hartwig conjecture 

We consider now the case when not only the phase S(t) but also the distribution function n(e) has multiple step-like 
singularities, Eq. (| A5|) . We will apply the generalized Szego formula Eq. (|A43[) derived above to this situation and 
then cut off the resulting logarithmic divergencies. While this is not a mathematically rigorous procedure, we should 
get in this way correct results for power-law exponents (apart from the summation over branches of the logarithms), 
cf. Sec. EH 

The relevant functions are now given by 



1 n 
3=0 



g®{z) = f >?S'rX.iUP*J 



(A44) 



with 



and 



N = ^- ( ln 



l + (e" 



l) aj 



2-rti \" 1 + (e~ iS ^ - l)a i+ i 
Here we have intriduced the notation 



U, 



-ln 



(i) 



3=0 



1 + (e-^ - l) aj 



1 + {e~ lS ' - l)oj- 



Li,j + ij+i j'+i — L 



(A45) 



(A46) 



2n 



log [1 + {e~ iSi - l)o 



The coefficients can be assigned to the singular points (ej,tj) in (e,t) plane, as illustrated in Fig. [6j 
Just as in Sec. I A II coefficients /3y control the asymptotic behavior of Uq at 5 — >• 00, 



3=0 



(A47) 



(A48) 
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Applying now Eq. (|A43[) with the cutoff procedure discussed in Sec I A H and representing the result in terms of physical 
variables, we find an asymptotic behavior of the determinant of interest 



A[<5(t),n(e)] = C[/3] exp 



i5Z*iAi(cj-A) 



0<i'<i<k 



j'<3 



A 



(A49) 



where we have used the condition 



(A50) 



i=0 



In full analogy with the case of Toeplitz determinants with Fisher-Hartwig singularities, see Sec. lA 11 Equation (IA49[) 
suffers from the ambiguity in the choice of /Jy . Pursuing this analogy, we conclude that the correct generalization of 
(|A49[) involves a summation over all possible choices of /3y , 



Affl,»i=E^ ex P 



br 



A 



■Etc 



0<i'<i<fe j'<j 

The summation over /3jj here can be understood as summation over all possible branches of logarithms in the definition 
of Lij, Eq. (|A47|) . with the additional constraints imposed [cf. Eq. (|A15p ] 



Loj — Lk+ij — Li <m +i = and Li t o — 



2vr 



A particular case of (|A51|) is the zero temperature determinant 



A[S{t),T = Q} = C [P] exp 



n [a(*» - *«')] 



(A52) 



(A53) 



0<i'<i<fc 



with ^ = (£; - <J i+ i)/27r. This result can also be checked by a direct calculation using the fact that at equilibrium 
the determinant is a gaussian functional of 5(t). 

Equation (|A51|) can be represented in a more symmetric way. Employng the condition 



_ Sj - S i+1 _ fl ( ) 
2tt 



(A54) 



(valid for arbitrary multi-step distribution) and combining it with (|A51|) and (|A53|) we arrive at the (cutoff indepen- 
dent) asymptotic of the normalized determinant A[5(t),n(e)] = A[5(t), n(e)]/A[6(t), T = 0], 



A[5(t)Xe)]=EC[/3]exp 

br 

where we have introduced the exponents 



■i ^2 Ufa 



i ,3 



0<i'<i<fe 0<j'<j<m 



Ifii' ,jj f PijPi* j' Pij 1 j 



(A55) 



(A56) 



Let us now discuss the applicability of the long time asymptotic expansion. Throughout the derivations of Eq. (|A55I) , 
it was assumed that the area of all the rectangles in the time-energy plane is large, i.e., for any fj, ij/, ej and 6j> 



{e l -e l ,)(t J -t-,) > 1, i^i', j^f. 



(A57) 



A straightforward extension of the cutoff procedure shows that asymptotics (|A55I) remains valid even if for some sets 
(i,i',j,f) the condition (IA57[) is not satisfied, provided one drops the corresponding factors from the product (|A55I) . 
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FIG. 7: The set of number /3y corresponding to the counting field <5(t), Eq. (|B1|) , in terms of j3?y and corresponding to 
individual pulses comprising 8(t). 



Appendix B: Factorization of the determinant for 5(t) consisting of several remote pulses 

In this section we discuss the determinant A[8(t), 11(e)} for the phase S(t) consisting of several rectangular pulses 
separated by large time intervals. We show that in this limit the determinant is given by a product of Toeplitz 
determinants corresponding to individual rectangular pulse. This question is relevant for the calculation of correlation 
functions in Luttinger liquid within generic boundaries (and in particular for "sharp boundary" model). In our current 
analysis we rely on Eq. (|A51| ). 

Let us consider counting field S(t) given by (see Fig. [Jj 

{0 = 5 , t < t 

Si, t < t < ti 

0, h<t<t 2 (Bl) 

S3, t 2 <t <t 3 

= S 4 , t 3 <t. 

Our aim is to show that as t 2 — ti goes to infinity, A[S(t),n(e) decouples into the product A[S^(t), n(e)] A[j( 2 ) (t), n(e)], 
with £W(t) and S^ 2 \t) being the two rectangular pulses constituting S(t). 

Let us introduce coefficients f)f) with I = 1,2 corresponding to the phases S (1 \t) and S^lt). Figure [7] demonstrates 

the connection between the coefficients /3.y corresponding to S(t) an d . Here we have introduced a set of integers rij 
satisfying rij = 0. Let us consider a particular term in the sum (|A51|) corresponding to the given set of coefficients 

ftfj and set of integers rij. The dependence of this term on t 2 — ti at t 2 — ti — > 00 is easily found to be 

(*a - t^^^M'^OS = (t a - t x )-^T=o^ . (B2) 
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FIG. 8: The real part of the determinant A[S(t),n(e)] for triple-step counting field S(t) (see inset) and triple-step distribution 
n(e). The solid blue line on both graphs represents the results of numerical evaluation of Fredholm determinant. Each curve 
corresponds to the fixed time ti (tz = 100, 200, . . . 500 for curves mark by letters from "a" to "e") while ti changes from to 
<2. Upper panel. Dashed lines show the fit of numerical result by expression (|C3() with coefficients Ci, C2 and C3 considered as 
fitting parameters. The same set of d is used for all curves. Lower panel. Dashed lines show the fit of numerical result with 
the correction terms (|C10[) taken into account. The oscillations of the determinant with time t\ are now correctly reproduced. 



We see that only terms of the sum (|A51[) characterized by rij = 0, for j = 0, . . . m survive the limit t2—t\ — > 00. Thus 
we find that in the long-time limit the determinant in question factorizes into a product of Toeplitz determinants 

A[5(t),n(e)}=A[6 {1 \t),n(e)}A[5 {2 \t),n(e)}, t 2 - h -> 00 . (B3) 

This factorization was proposed in Ref. 19,20 on the basis of physical arguments and employed to calculate single- 
particle Green functions. In the main text of the present paper we also use an analogous factorization for the case of 
two-particle correlation functions. 



Appendix C: Numerical check for the asymptotics (|A55|) of determinants with multiple t and e discontinuities 



In this section we present a numerical verification of the asymptotic formula (|A55[) for determinants with multiple 
discontinuities both in A{t) and in B(e). We calculate numerically the normalized Fredholm determinant A[S(t),n(e)] 
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with a triple-step counting field (see inset of FigJH]) 

8(t) = 



0, 

51 = 1.1 x 2tt, 

5 2 = 0.98 x 2tt . 
0, 



f < 

o < t < h 
ti <t <t 2 
t 2 <t, 



(Cl) 



as a function of t\ and t 2 and compare results to the predictions of Eq. (|A55[) . The energy distribution is taken to be 



n(e) 



e < e 
e < e < ei 
ei < e < €2 
£2 < e , 



(C2) 



with eo = —3/4, ex = —1/2 and £2 = 1/4. The numerical procedure used here is based on the time discretization and 
is analogous to the one employed in Ref. 23 in course of studies of Toeplitz determinants. 

Taking into account that Si and S 2 in (ICip are close to 2-7T, one infers the following three dominant contributions 
to the sum (|A55|) : 



-d) 



(C ie - ie °* 2 + C 2 e-" 2t2 ) tl 1 f 2 2 (t a - h)^ + Cae-^HJ 1 f 2 2 (t 2 - h)^ 



Here we have absorbed the dependence of the determinant on energies €j into the coefficients 
The exponents in (|C3[) are given by 



-L(S 1 -2w)(S 1 -5 2 ) 
2tt z 



7 2 



(i) _ 



-^(-£i<S 2 + 2tt(5i + 2vr(52 - 4tt 2 ) 

27T Z 



712 



^(<5 2 -2^)(<5 1 -5 2 ) 



1 

2~ 
1 

2~ 



27T 

72 2) = ttt(-^ S 2 + 4tt5i + 47r<5 2 - 12tt 2 ) 

Z7T 



7f 2 } = ^ (<S 2 - 4^(^-52). 

Z7T 



(C3) 

(C4) 
(C5) 
(C6) 
(C7) 
(C8) 
(C9) 



A characteristic feature of the exponents (|C4I) — (|C9[) is their smallness at <5i ~ ~ 27r. 

The upper panel of FigfSJshows the real part of the determinant A[S(t),n(e)}. The solid blue lines on both graphs 
represent the results of numerical evaluation of Fredholm determinant. Each curve corresponds to the fixed time 
t 2 (t 2 = 100, 200, ... 500 for curves mark by letters from "a" to "e", respectively) while t\ changes from to t 2 . 
Dashed lines show the fit of numerical result by expression (IC3[) with coefficients C% , C 2 and C3 considered as fitting 
parameters. The same set of Cj is used for all curves. We see that the fit reproduces correctly the overall behavior of 
the determinant but the oscillations with time t\ are not captured in this approximation. An analysis of (|A55|) shows 
that for chosen parameters the dominant correction to (IC3[) is given by 



A 



<2) 



(a 



with 



4 e 



(3) 

7i 



-ie n *2— j£2*i 



c. 



-ie 2 t 1 -ie 1 t 2 



(3) (3) (3 > 



2n 

(3) 
7 2 



MS2-S 1 ) + -(35 1 - 25 2 )-3.. 



2tt ; 



r ^ 2 + -(25 2 + ^)-3, 



7l2 



2tt 2 



6^ 



(C10) 

(Gil) 
(C12) 
(C13) 



Taking the correction (jClOl) into account, we obtain a fit to numerical data (bottom panel of Fig. [S]) which correctly 
captures the oscillation in time. The resulting agreement is essentially perfect. We thus conclude that our conjecture, 
Eq. (IA55[) . is fully supported by the numerical simulations. 
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